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Existence and order of the phase transition 
of the Ising model with fixed magnetization 
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Properties of the two dimensional Ising model with fixed magnetization 
are deduced from known exact results on the two dimensional Ising model. 
The existence of a continuous phase transition is shown for arbitrary val- 
ues of the fixed magnetization when crossing the boundary of the coexis- 
tence region. Modifications of this result for systems of spatial dimension 
greater than two are discussed. 
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1. INTRODUCTION 

In a recent paper, 2d and 3d Ising models with fixed magnetization are 
investigated numerically and signatures are found in the microcanonical 
caloric curves of finite systems which hint a first-order phase transition. 
From the numerical data, however, it is not clear if the observed signa- 
tures scale appropriately to persist in the thermodynamic limit, and hence 
properties of the infinite system cannot be inferred from the data. 

Although an exact solution is known only for the zero-field case of 
the 2d Ising model without any constraints on the magnetization,^ there 
are a number of exact results which can be used to tackle the questions of 
existence and order of a phase transition at fixed magnetization analyti- 
cally. As a convenient thermodynamic function to discuss this topic, the 
entropy as a function of the interaction energy and the magnetization is 
chosen. 

An important ingredient in the following discussion is the equiv- 
alence of ensembles, which holds for Ising systems of arbitrary spatial 
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dimension.® This allows to combine the known exact results, which are 
typically obtained in the canonical ensemble, with some simple geometri- 
cal arguments on the microcanonical entropy. 

Sec. H gives a definition of the Ising model. In Sec. |3|, several ther- 
modynamic quantities are defined and implications of some known exact 
results of the 2d Ising model on these quantities are discussed. Then, only 
elementary analysis is needed to establish the existence of a phase transi- 
tion of the Ising model with fixed magnetization in Sec. f| and to identify 
the transition as a continuous one in Sec. |5|. Modifications of this results 
for the case of spatial dimension d > 2 are discussed in Sec. ^ 

2. THE 2D ISING MODEL 

Consider an evenQ number N = L 2 , L G 2Z, of classical spins 
Oi G {— 1,+1}, i = 1, ... ,JV, on a two dimensional quadratic lat- 
tice {1,... , L} 2 with periodic boundary conditions. Then the nearest- 
neighbor Ising modem is defined by the Hamiltonian 

H:T N ^R, a (-»• £{a) - hM{a) (2.1) 

where = {— 1,+1} is the configuration space of the system, a = 
(ax,... ,0~n) £ Tat are called configurations, and h G R is an external 
magnetic field. 

N 

M : T N -»• (2Z) n [-N,+N], a^^cTi (2.2) 

is the magnetization and 
£ : T N -> (4Z n [-2N, +2N]) \ {-2N + 8, 2N - 8}, a ^ 

(id) 

(2.3) 

is the interaction energy. denotes a summation over all pairs of spins 
which are neighbors on the lattice. 

3. THERMODYNAMIC QUANTITIES 

We follow Ref. ® to define the entropy density^ s of the 2c? Ising model in 
the thermodynamic limit iV — >• oo as a function of the interaction energy 

2 This is only for notational simplicity. 

3 The term "density" will be omitted in the following 
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density 2 e and the magnetization density 2 m. Let 

B e (e,m) := {(e',m') G R 2 : (e - e 1 ) 2 + (m - m f < q 2 } (3.1) 
be the ball of center (e, m) and radius Then the entropy can be defined 



as 



s : V{s) -> R, (3.2) 



min lim — In 

0>O JV-i-oo TV 



where | • | denotes the cardinality of a set. The domain 

V{s) := |(e,m)|e G [-2, +2], |m| < (3.3) 

of the entropy has the shape of a triangle. This is due to the fact that 
a given value of the interaction energy e is equivalent to a given propor- 
tion of "antiparallel" neighboring spins Oi = —<jj neighbors), which 
implies an upper bound on the absolute value of the magnetization for 
configurations of the given energy. Due to the spin inversion symmetry 
of the Ising Hamiltonian, the entropy s(e,m) is symmetric with respect 
to the magnetization m. Fig. ffl is intended to give an idea how the graph 




Figure 1: Approximate graph of the entropy as a function of the interaction 
energy e and the magnetization m. 
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of s approximately looks like. An exact closed form expression for s is 
not known. Via Legendre transformation, this would be equivalent to a 
closed form expression of the free energy density of the 2d Ising model for 
arbitrary external field h. 

The variable thermodynamically conjugate to the magnetization m 
is /3h, the product of the inverse temperature (3 and the external field h. 
An expression of the magnetization 

m: [-2, +2] x (R\{0» -»• [-1,+1], (e, /3h) i-> m(e, (3h) (3.4) 

as a function of f3h and the interaction energy e is obtained from the 
entropy implicitly via the Legendre transformation 

sup[s(e, m) + f3hm] =: s(e, fh(e, (3h)) + /3hm(e, (3h) (3.5) 

m 

The spontaneous magnetization m± is defined^ as the zero-field limit 

m± : [-2, +2] -> [-1,+1], £^ lim m{e,(3h) (3.6) 

/3/i— >±0 

The domain of the entropy s, the triangle will be used in the 

schematic picture of Fig. |^ to illustrate some properties of the entropy. 
The bold curve in the triangle is the spontaneous magnetization of the 
2d Ising model, for which an exact expression is known from combining 
Yang's resu with the caloric curve obtained from Onsager's solution.™ 
For the spontaneous magnetization the condition 

m± (e) = ^ e > -V2 (3.7) 

holds, and the value —a/2 is called the critical energy. 

The region of T>(s) bounded by the non-zero spontaneous magnetiza- 
tion and the line e = — 2, 

C:=|(£,m)GD(s) \m\ < |m±(e)|J- (3.8) 

is called the coexistence region (hatched in Fig. ||). Inside the closure C 
of the coexistence region, s is constant with respect to the magnetization, 
i.e., 

s(s,m) = s(e,m±(e)) W(e,m) eC (3.9) 



4 This is not entirely in accordance with the standard terminology, where one speaks 
of spontaneous magnetization only where m±(e) ^ 0. 
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Figure 2: Schematic picture of the entropy s(e, m). The triangle is the domain 
of the entropy, the bold line is the spontaneous magnetization. Parallels drawn 
inside the coexistence region C (hatched) indicate straight lines of constant 
entropy. 



which is the signature of the field-driven first-order phase transition of the 
2d Ising model at low energies and zero external field. In Fig. [^, some 
parallel lines are drawn inside the coexistence region to symbolize lines of 
constant entropy. Outside the coexistence region C, the concavity of the 
entropy combined with the definition of the spontaneous magnetization 
fl3.6|) implies the strict inequality 



s(e,m) < s(s, m±(e)) V(e,m) ^ C, m ^ 



(3.10) 



Now, from the two-variable function s, we define the one-variable function 



s zf : £>(s zf ) 



R, 



e I— > s[s, m: 



(3.11) 



m=m± (e) 

on the domain T>(s z {) = [—2, +2], which is the entropy of the 2c? Ising 
model in the zero-field limit, and the family of one- variable functions 



R, 



171=1710 



(3.12) 



on the domain V(s mo ) = [—2,2 — 4|mo|], which is the entropy of the 2d 
Ising model with fixed magnetization itlq. In the following, existence and 
order of a phase transition in the 2d Ising model with fixed magnetization 
are discussed by making use of these two functions. 



6 



Michael Kastner 



4. EXISTENCE OF A PHASE TRANSITION 

For all values of the fixed magnetization |m | < 1, the entropy s mo is 
non-analytic for at least one value of e G T>(s mo ). This is a consequence 
of: 

a) The fact that every line {(e, m) G int'D(s) \ m = m } of constant 
magnetization itlq G (— 1,+1) in the interior of T>(s) has at least 
one point of intersection with the spontaneous magnetization (as 
sketched in Fig. 0), 



b) s TOo (e) = Szf(e) for all (e, mo) G C [from Eq. (J379f) l , 

c) s mQ (e) < s zi (e) for all (e, m ) ^ C, m 7^ [from Eq. ( [3.10| )1, 



d) The entropy s z f(e) of the Ising model in zero field is analytic for all 



ds zt (e) 
de 



e < —y/2 and non-analytic for e = — y/2 [Otherwise (3(e) 
would show a non-analyticity for e < — y/2. Due to the equivalence 
of ensembles, this would be in contradiction to the analytic result 
for the caloric curve s{(3\ below the critical temperature obtained 
from Onsager's solution.®], 

e) s mo=0 (e) = s z f(e) for all e [from Eqs. ( |3.7| ) and (|3.9|)1. 

Case m 7^ 0: Item (a) guarantees that, for all values of the fixed 
magnetization |mo| < 1, the line of constant mo crosses the boundary 
of the coexistence region once, and the corresponding energy £o(mo) is 
given by the inverse function of the spontaneous magnetization m±. From 
conditions (b) and (c), we know that s z { and s mo are identical on the 
interval [— 2, £ (mo)], but different elsewhere. According to (d), s z f is 
analytic on the interval (—2, —y/2), and, as the analytic continuation of a 
function is unique, it follows that s mo (e) is non-analytic at e = Eo(m ) < 
—y/2 for all < |m | < 1. 

Case m = 0: In this case, according to item (e), the entropy func- 
tion of the Ising model with fixed magnetization is identical to that of 
the zero field case on the entire interval [—2, +2]. Hence, also the non- 
analyticity of s z f(e) at e = e (0) = —y/2 [item (d)] is present in s mo=0 . 

From Yang's result® for the spontaneous magnetization of the 2c? 
Ising model, the inverse temperature (3q corresponding to £o(mo) via the 
caloric curve is known to be 

fa : [-1,+1] -»• [/? c ,oo], m i-> iarcsinh (l - m s )~ 1/A (4.1) 
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where (3 C = | ln(l + \fi) is the critical temperature of the 2d Ising model. 
Hence, we conclude that the 2d Ising model with fixed magnetization mo 
shows a phase transition at the inverse temperature /?o(^o)- 

5. ORDER OF THE PHASE TRANSITION 

In this section, we show that the phase transition of the Ising model with 
fixed magnetization is a continuous one. 

A temperature-driven first-order phase transition is characterized by 
the appearance of latent heat, i.e., a discontinuity in the caloric curve 
e(f3), where /3 is the inverse temperature. In the Ising model with fixed 
magnetization, due to the thermal equation of state 

/3(E , ro) = ^) (5.1) 

such a discontinuity corresponds to a linear piece in the entropy s mo . In 
the following, by proving the absence of a linear piece in s mo , the phase 
transition of the 2d Ising model with fixed magnetization m is shown 
to be continuous for all values of m G [— 1,+1]. For this discussion we 
distinguish between the two cases of a linear piece inside and outside the 
coexistence region, respectively. 

5.1. No linear piece of s mo inside the coexistence region 

From Eq. ( |3.9| ) it follows that a linear piece of s mo inside the coexistence 
region implies a linear piece in s z {. It is known from Onsager's solution® 
that there is no temperature-driven first-order phase transition in the 2d 
Ising model in zero field, hence there is no linear piece of s mo inside the 
coexistence region. 

5.2. No linear piece of s mo outside the coexistence region 

This result can already be found in the literature. For an Ising system of 
arbitrary spatial dimension in non-zero external field or at temperatures 
above the critical temperature, there exists a unique, translation invariant 
equilibrium state.® It is shown in Ref. ® that this implies strict concavity 
of the entropy outside the phase coexistence region C. This strict con- 
cavity of course does not allow for a linear piece in the entropy s mo of an 
Ising system with fixed magnetization m outside C. 

Additionally, another proof of the absence of a linear piece in s mo 
outside the coexistence region is given by showing a stronger result on the 
analyticity properties of the entropy s. It will be argued that a linear piece 
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of s mo outside the coexistence region is in contradiction to the analyticity 
properties of the Gibbs free energy 

g : V(g) -> R, (/?, (3h) ^ - lim {(3N)- 1 In V (5.2) 

TV— >oo ' J 

T>(g) := R + x R, of the 2d Ising model. It follows from the circle theorem 
of Lee and Yang® and a result from Lebowitz and Penrose®* on analyt- 
icity properties, that g is analytic for all non-zero values of the external 
field h or for all inverse temperatures (3 smaller than the critical inverse 
temperature.^ The Gibbs free energy g and the entropy s are connected 
by means of a Legendre transformation 

s(e, m) = sup [—(3g{(3, f3h) + (3e — (3hm] (5.3) 

f3,/3h 

For regions where g is differentiate, e and m can be expressed in terms 
of the Gibbs free energy as 

and 

e((3, fih) = g((3, (3h) + p dg ^ h) (5.5) 

For regions where g is strictly concave, the inversions (5{e, m) and (3h(e, m) 
exist, and Eq. (|5-3) can be rewritten as 

dg(J3(e,m),ph(e,m)) 



s(e, m) = (3(e, m) 



(3(e, m 



+(3h(e, m) 



d(3(e, m) 
dg(J3(e,m),/3h(e,m)) 



(5.6) 



d/3h(e, m) 

This implies analyticity of s(e,m) for all pairs (e(f3, [3h),m(f3, (3h)) G P 
for which g(/3,/3h) is analytic, and 

P = e V(s)\c\ > 0} (5.7) 

is the set of pairs (e, m) outside the closure of the coexistence region 
corresponding to positive temperatures. Then, also s mo (e) is analytic for 
all pairs (e, m ) G P. 

The entropy s is zero on the boundary of its domain and continu- 
ously differentiable on its entire domain^ . Together with the analyticity 
properties discussed above, this rules out the possibility of a linear piece 
in s mo outside the coexistence region. 

5 Actually, this result is valid for arbitrary spatial dimension. 
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6. SPATIAL DIMENSION GREATER THAN TWO 

Some remarks are in order concerning the validity of the arguments of 
Sees. |] and |5| for Ising models of spatial dimension d > 2. The crucial 
difference to the 2d Ising model is that, to the knowledge of the author, 
it is not clear if there exists more than one temperature-driven phase 
transition in the Ising model in zero field for d > 2. If this is not the 
case, the results of this paper are valid for Ising systems of higher spatial 
dimension. If there is more than one such transition, some modifications 
of the results on the Ising model with fixed magnetization arise: 

• For values of (e, m ) inside the coexistence region C, the Ising model 
with fixed magnetization shows identical thermal behavior (in the 
sense of identical caloric curves) as the Ising model in zero field, 
including phase transitions possibly existing in the zero-field case. 



The results used in Sec. are valid for Ising systems of arbitrary 
spatial dimension. Hence, there is no linear piece in the entropy 
of the Ising model with fixed magnetization outside the coexistence 
region C . 

In general the arguments of Sec. f| are also valid for d > 2 and 
the Ising model with fixed magnetization shows a phase transition 
when crossing the boundary of C . However, when crossing this 
boundary at an energy for which the zero-field entropy s z { is non- 
analytic, it might be the case that the entropy of the Ising model 
with fixed magnetization is the analytic continuation of s z { across the 
locus of the non-analyticity. Then, for this particular value of the 
fixed magnetization mo, the Ising model with fixed magnetization 
does not show a phase transition when crossing the boundary of the 
coexistence region. 



7. CONCLUSION 

In the 2d Ising model with fixed magnetization, for all values of the 
fixed magnetization |m | < 1, there is a continuous phase transition. 
This phase transition is unique in the sense that there are no further 
non-analyticities in any thermodynamic function or in the caloric curve. 
The inverse temperature at which the phase transition occurs is given by 

A)(^o) — \ arcsinh (1 — wig) ^ 4 . Hence, in the thermodynamic limit, 
the S-shape reported for the microcanonical caloric curve of finite Ising 
systems with fixed magnetizations in Ref. 1, does not converge to a straight 
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line signalling a first order phase transition, but to a single point at which 
a continuous phase transition occurs. 

For spatial dimensions d > 2, a first-order transition can occur in the 
Ising model with fixed magnetization only if there is a first order transition 
in the <i-dimensional zero-field case. Although, to the knowledge of the 
author, a proof is lacking, this is expected to be not the case for the 3d 
Ising model. 
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